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ABSTRACT 

The First Ionization Potential (FIP) effect is the by now well known enhance- 
ment in abundance over photospheric values of Fe and other elements with first 
ionization potential below about 10 eV observed in the solar corona and slow 
speed solar wind. In our model, this fractionation is achieved by means of the 
ponderomotive force, arising as Alfven waves propagate through or reflect from 
steep density gradients in the solar chromosphere. This is also the region where 
low FIP elements are ionized, and high FIP elements are largely neutral leading 
to the fractionation as ions interact with the waves but neutrals do not. Helium, 
the element with the highest FIP and consequently the last to remain neutral 
as one moves upwards can be depleted in such models. Here, we investigate this 
depletion for varying loop lengths and magnetic field strengths. 

Variations in this depletion arise as the concentration of the ponderomotive 
force at the top of the chromosphere varies in response to Alfven wave frequency 
with respect to the resonant frequency of the overlying coronal loop, the magnetic 
field, and possibly also the loop length. We find that stronger depletions of He 
are obtained for weaker magnetic field, at frequencies close to or just above the 
loop resonance. These results may have relevance to observed variations of the 
slow wind solar He abundance with wind speed, with slower slow speed solar 
wind having a stronger depletion of He. 

Subject headings: Sun:abundances - Sun: chromosphere - turbulence - waves 



Introduction 



Our k nowle dge of the solar wind has come a long way since the original prediction 
of iParkerl (119581 ) . It is now known to have two separate components; "fast wind" which 
originates in coronal holes with speeds generally greater than 600 km s _1 , and "slow wind" 
with speeds less than 600 km s _1 , the precise solar origin of which is uncertain. The fast 
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solar wind is usually rather stead y in speed and mass fl ux, and is generally taken to be the 
solar wind of Parker's prediction (lAntiochos et al.ll201l[ ). although it is now rec ognized to be 
accelerated close to the sun by the damping of ion cyclotro n and MHD waves (jOfmanll2010l ) 
rather than by heat conduction as in the original concept (|Parkerlll963l ). 



By contrast the slow solar wind is unsteady, with variability presumably reflecting solar 
source regions subject to varying phenomena associated with magnetic activity on the Sun. 
The slow solar wind is also distinguished from the fast by the properties of its turbulence. 
The fast wind exhibits essentially purely Alf venic turbulence, wherea s the slow wind also has 



magnetic structures advected within it (e.g. iBruno fe Carbondl2005l ). Since such structures 



may also be represented by superpositions of counter-propagating Alfven waves in equal am- 
plitude, turbulence in the slow wind is sometimes said to be "more balanced". The elemental 
composition of the two types of solar wind is also known to be different. Solar plasma in 
closed coronal loops and in the slow speed wind is known to be subject to a fractionation 
process, as the material that supplies these features moves up from the chromosphere to 
the corona. Elements which are strongly ionized in the chromosphere like Fe, Si, Mg, with 
first ionization potentials less than about 10 eV (low FIP elements), are enhanced in abun- 
dance by a factor of about 3-4. High FIP elements are generally unaffected. He often shows 



a depletion in the solar corona ^Laming fe Feldman 2001 . 2003) and wind , which is larger 



and more variable in t he slow wind ( lAellig et al. 



2001 



Kasper et al.l 120071 ) than in the fast 



( IReisenfeld et al.ll200ll ). The abundance enhancement of low FIP ions in the fast solar wind 
is also much smaller than in the closed loop corona or the slow wind. 



This variability arises naturally in models of the FIP effect (see iLamingl 120041 . 12009 



20121 . and discussed in more detail below), where the agent of fractionation is the pondero- 
motive force associated with Alfven waves propagating between the chromosphere and the 
corona. The changing magnetic field geometry between open and closed field lines alters 
the propagation of Alfven waves in such a way that the FIP fractionation produced by the 
ponderomotive force varies in the manner observed, i.e. more fractionation on closed field 
lines. Thus we assume henceforward that one fundamental difference between fast and slow 
solar wind is that slow solar wind must have originated on a closed coronal loop to become 
fractionated, and that this loop is subsequently opened up to release its plasma into wind. 

Some authors argue that this distinction is not required. While the slow wind speed is 
well k nown to correlate inver sely with the expansion rate of open field in equatorial coronal 
holes (jWang fe Sheeleylll990l ) leading to the association o f the slow wind with super-radially 
expanding flux tubes at the edges of coronal holes (e.g. iKo et al.ll2006l ). th e view that the 
slow w ind originates exclusively in open field regions has been advanced by ICranmer et al. 
(120071 ) . They use their simulated Alfven wave spectra in polar and equatorial coronal holes to 
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calculate FIP fractionation following iLamingl (120041 ) and find a difference between fast wind 
from polar holes and slow wind from equatorial hole similar to that that is observed. On 
this basis, the transition in element al composition from fast to slow wind sho uld probably be 
gradual. However in observations (IGeiss et al.l Il995l ; IZurbuchen et al.l 119991 ) this transition 
is abrupt, both in th e change in compositio n and the freeze- in temperature associated with 
the O charge states. lAntiochos et al.l ( 120111 ) also comment that a tight correlation between 
element abundances and wind speed should result from such models, a finding which also 
is not supported by observation. We therefore adopt the view that the fast and slow winds 
are distinct entities, and the distinction arises as the fast wind originates exclusively in open 
field, while to become fractionated, the plasma comprising the slow wind must have spent 
time inside a closed coronal magnetic loop, which is subsequently opened up by reconnection. 

In th is paper we wi ll investigate the abundance pattern expected from such a process. 
We follow lLamingi ( 120121 ) in assuming the Alfven waves to have a coronal origin, presumably 
associated with heating mechanisms such as n anoflares and their associate d reconnection (e.g. 
Sturrockill999l ) or resonant abso rption (e.g. iRuderman fc Roberts! 120021 ) . Due to repeated 
reflection at each loop footpoint ( jOfmanl 120021 ). equal intensity Alfven waves propagating in 
each direction are expected to be produced by these processes. We pay particular attention 
to variations of the abundance of hel ium, observed to change with slow solar wind speed 
( Aellig et al. 2001 ; Kasper et al. 2007|). We show that not o nly does the FIP effect model 
based on the ponderomotive force (lLamingi 120041 . 120091 . 120121 ) explain the fact that the solar 
wind abundance of helium is lower than its generally accepted photospheric value, but that 
variations in this abundance may also be within reach of the model. While processes such 
as gravitational settling or inefficient Coulomb drag are not completely ruled out, their 
importance in forming the elemental composition of the solar wind is reduced from that 
discussed in previous works (IBochslerl 120071 ; iKasper et al.l 120071 ). 



2. The Ponderomotive Force Model 



In the model developed by lLamingi (120041 . 120091 . 120121 ). FIP fractionation has its origin 
in the action of the ponderomotive force on chromospheric ions, but not neutrals. The 
ponderomotive force is a second order term incorporating the p5v ■ V5v and 53 x 5B/c 
terms in the MHD momentum equation. In many cases it looks like a force due to the 
gradient of wave magnetic pressure, but this obscures its true nature as the interaction of 
waves and particles through the refractive index of the plasma. Given this, it can be seen 
that both upward and downward ponderomotive accelerations may occur, with the latter 
giving a potential explanation for the so-called "inverse FIP" effect seen in active stars and 
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in dwarfs of later spectral types than the sun (e.g. IWood Linskyil2010l ). The fractionation 
arises from the time averaged part of the ponderomotive force, which for waves of frequency 
ua << fi, the ion gyrofrequency, is given by (e.g Litwin Sz Rosner 1998, Laming 2009) 
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where q is the ion charge. The ponderomotive acceleration, F/m, is independent of the ion 
mass. 



Our basic approach builds on iHollwed (119841 ) , and is illustrated by a cartoon in Figure 
1. We start with a downward propagating Alfven wave at the /3 = 1.2 layer of one footpoint 
of a coronal loop, and integrate the non-WKB transport equations spatially for the forward 
and backward propagating Alfven waves through the corona to the other footpoint, where 
we now have a mixture of waves propagating in both directions, depending on the resonant 
properties of the loop with respec t to the wave frequency. The transport equations are 
(jCranmer &: van Ballegooijenll2005l . and refs therein) 
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where z± = 5v±±5B±/ ^/A-np are the Elsasser variables for Alfven waves propagating against 
or along the magnetic field respectively. The Alfven speed is Va = B/y/Airp, the upward flow 
speed is u, the mass density is p, and the magnetic field is B. The signed scale heights are 
Hd = p/ (dp/dr) and Ha = VaI {dV^/dr). In the solar chromosphere and corona u « Va, 
and we put u — 0. Once z± are determined by numerical integration of equation 2, we 
calculate 5v± and SB^/y/Airp, and hence the perpendicular electric field SE±_ from 
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Contributions from waves of the same frequency propagating in both directions along the 
magnetic field are naturally accounted for. In this work, for reasons of pedagogy, we only 
consider one pair of parallel forward and backward propagating Alfven waves of a given 
frequency at a time, and neglect any nonlinear effects or dissipation. We are tacit ly assuming 



that l oop dynamics are mainly driven Alfven wave interactions, as suggested by iBigot et al. 
j2008h . 



The degree of fractionation is given by the formula 

fractionation = exp ^2 J ^ s a{y e fjjv s ^i)lv 2 s dz 



(4) 
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(see equation 9, Laming 2009, equation 12, Laming 2004, which follow Schwadron et al. 
1999), where £ s is the ionization fraction of element s, a is the ponderomotive acceleration, 
u eff = v s,i v s,nl \^s v s,n + (1 — £,s) ^s,i] where v s i and z/ Sj „ are the collision rates of ions and 
neutrals, respectively, of element s with the ambient gas. Also v 2 s = kT/m s + v 2 turb + v 2 urb , 
where v^turb is the amplitude of microturbulence coming from the chromospheric model, and 
v turb is the amplitude of longitudinal waves induced by the Alfven waves themselves. These 
derive from the harmonic variation of the ponderomotive force, and unlike the fracti onation, 
do no t require an external density gradient. A full description of this is given in lLaming 
(120121 ). The limits of integration, z\ and z u bracket the region where the ponderomotive 
acceleration is effective in the chromosphere. Typically z\ is at the /5 = 1.2 layer in the low 
chromosphere, and z u is in the transition region where all elements are ionized. 



The chromosphere at each footpoint can be based on any o f thelVernazza. Avrett. fc Loeser 

( 1981 ) models or similar. Here we use the Avrett fc Loeser ( 20081 ) update of VALC. The 
ionization balance of the minor ions is computed at each height in the chromosphere using 
the model temperature and electron density, and a co ronal UV- X-ra y spectrum appropri- 
ately absorbed in the intervening chromospheric layers. iLamingl (120121 ) investigated different 
approximations for computing the chromospheric ionization balance, and while the degree 
of FIP fractionation varied somewhat, the validity of this explanation for the phenomenon 
remained intact. The chromospheric magnetic field (illustrated in Figure 2) is taken to be a 
2D force free field from Athay (1981) and designed to match chromospheric magnetic fields 
in Gary (2001) and Campos & Mendes (1995), which represents the expansion of the field 
from the high /3 photosphere where the field is concentrated into small network segments, 
into the low (3 chromosphere where the field expands to fill much more of the volume. The 
coronal section of the loop has constant magnetic field, and a density scale height of 125,000 
km. 

The precise model for the chromospheric magnetic field is not critical, because most of 
the FIP fractionation occurs in the upper chromosphere where the field is most uniform and 
vertical. In the lower chromospheric region with curving magnetic field in the plasma layer 
where sound and Alfven speeds converge, contrary to the case in the corona, wave interaction 
processes play an important role. We do not attempt to model these directly, but argue that 
the increased longitudinal wave amplitudes in this region, arising from mode and parametric 
conversion, increase the value of v 2 in the denominator of the integrand i n equati o n 4 t o 
such an extent so as to render fractionation insignificant in this region. In ILamingl (120121 ). 
the longitudinal and Alfven wave inten sities are coupl e d by equat ion 15 in that ref erence, 
motivated by calculations il lustrated inlCranmer et al.l (120071 ) and iKhomenkol (120101 ). Here 
we note that recent work by lHeggland et al.l (120111 ) suggests that sound waves may be strong 
in this region even without the presence of Alfven waves, as a result of turbulent motion in 
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the convection zone. 



3. A Parameter Survey of FIP Fractionation 



In this paper we begin to explore the parameter space of magnetic field and loop length 
looking at coronal mag netic fields of 5 G, 10G, 15G, and loop lengths of 50,000 km, 75,000 
km, and 100,000 km. ILamingi ( 120121 ) investigated the variation of FIP fractionation in 
a 100,000 km loop with 20G magnetic field, giving a loop resonance angular frequency 
of 0.07 rad s _1 , as the Alfven wave driver of the FIP effect varied with respect to this 
resonance. Wit h wave amplitudes c hosen to give a fractionation Fe/O of about 4 as observed 
for example by lBryans et al.l (120091 ) . on resonance, He/O was about 0.6, S/O about 1.7, and 
C/O about 1.2. As the wave frequency increased above the resonant value, as might happen 
for example if a coronal heating process excited waves at the loop resonance, followed by 
chromospheric evaporation reducing this resonant frequency, He/O decreased slowly at first 
and then increased up towards photospheric values, S/O increased rapidly then levels off 
to a maximum of ~ 1.9, and C/O increased rapidly to a maximum of 1.6 and then began 
decreasing. Other fractionations varied also, but these three appeared the most robust and 
easiest to interpret in terms of available observational data. 

Performing similar calculations for a wider variety of loop models, again keeping the 
Fe/ O fractionation to a value of about 4, we see a similar pattern with Alfven wave frequency 
as was seen for 20G and 100,000 km. In Figure 3 we plot an example of how the fractionation 
varies with Alfven wave frequency, for a field of 15 G and a loop length o f 75,0 00 km. Overlaid 
on the plot are the fractionations measured in the Ivon Steiger et al.l (120001 ) Ulysses survey 
of the solar wind. The color bar gives the range of frequencies examined, which start from 
just below resonance at 0.651 rad s -1 . The black line is below resonance. The purple line 
is the resonant case. The two blue lines that still have strong helium depletion are just 
above resonance. The green, yellow and red lines are further from resonance, have less and 
less helium depletion and enhanced S and C. In general the modeled fractionation pattern 
between Mg, Fe, and Si does not match the observations, with Fe predicted to have more 
FIP fractiona tion than actually observed in the solar wind. However oth er spectroscop ic 
observations (IBryans et al.l 120091 ) actually match this pattern very well (see ILamingi 120121 ). 



In Figure 4, we plot the He/O abundance ratio as a function of frequency, in units of the 
loop resonant frequency, for a variety of loop lengths and magnetic field strengths. The curve 
for 100,00 km loop length and 20 G magnetic field is very similar to ILamingi (120121 ). with 
a minimum He/O fractionation of 0.6 at resonance. The position of this minimum shows 
a clear trend to lower He/O and higher frequency above resonance for weaker magnetic 
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fields (and possibly also longer loop lengths); cases where the ponderomotive force is more 
concentrated at the top of the chromosphere. Figure 5 gives the coronal amplitude of the 
Alfven wave in each case. For 5G and 100,000 km, and 10G with loop lengths 75,000 and 
100,000 km it was not always possible to achieve a fractionation for Fe/O of about 4 for any 
coronal wave amplitude, so these points were omitted in Figures 4 and 5 and in Figures 9 
and 10 to follow below. The coronal wave amplitudes at resonance for each case are very 
similar, suggesting that if we normalized by this parameter instead of forcing Fe/O ~ 4, we 
would have very similar results. The curve for 15 G, 75,000 km (the upper green curve) gives 
the coronal wave amplitudes for the fractionations shown in Figure 2. 

Figure 6 shows the energy transmission characteristics of the various loops considered 
here, calculated from the fluxes of upgoing and downgoing waves in the loop footpoint 
encountered at the end of the numerical integration, and displayed with the same color 
coding as in preceding figures. On resonance, all loops have an energy transmission of unity. 
The width of this resonance increases with decreasing loop length, as does the minimum 
value of the energy transmission in between resonances. A trend with loop magnetic field is 
less easy to discern. 

Figure 7 depicts the main features of the simulation for the case of 15 G, 75,000 km and 
0.07 rad s~ 1 just above the resonant frequency of 0.0651 rad -1 . Top left gives the density 
and temperature structure of the chromosphere to be read on the left (black) and right hand 
(gray) axes respectively. Bottom left gives the ponderomotive acceleration and slow mode 
wave amplitude. The "spike" in the ponderomotive acceleration can be seen to coincide 
with the steep density gradient at the top of the chromosphere. The ionization structure in 
the chromosphere, illustrated in the top right panel for Fe, Mg, S, and C, and also in the 
bottom right panel (gray lines, to be read on the right hand side axis) along with the FIP 
fractionations (black lines, to be read on the left axis), determines which elements can be 
accelerated and therefore fractionated. 

Fe, Mg, S and C are all mostly ionized throughout the relevant region of the chromo- 
sphere, H and O become ionized between 1500 and 2000 km altitude in the chromosphere, 
and He is only ionized at the very top of the chromosphere. With our normalization of Fe/O 
~ 4, He becomes depleted relative to O when the ponderomotive acceleration is concentrated 
at the top of the chromosphere, i.e. when the "spike" associated with the steep density gra- 
dient is very strong, as in the case illustrated. O is accelerated up into the chromosphere 
while He is left behind, being the last element in the chromosphere to remain neutral. When 
waves penetrate deeper in the chromosphere Fe, Mg, S and C are all easily fractionated, and 
the relative contribution from the top of the chromosphere where O is accelerated and He 
left behind is less important. Since we are modeling specifically for Fe/O of 4, He/O will 
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be lowest when the acceleration region is confined to the top of the chromosphere, which 
happens close to the loop resonant frequency. 

There are several factors that play into how deep waves penetrate into the chromosphere. 
Longer wavelength waves resonant with longer loops are not expected to penetrate so far into 
the chromospheric density gradients before being reflected back into the loop. This is clear 
from the fact that for a given loop, the first harmonic frequency is slightly less than twice 
that of th e fundamental. Similar results are known for photospheri c Alfven waves launched 



upwards (ICranmer fc van Ballegooijenl 120051 ; IVerdini fc Veilil 120071 ) . However in our case 



this effect is probably rather small, and not discernable on Figure 4. Figure 8 illustrates 
this in more detail. The four panels show at top left the real (black) and imaginary (gray) 
parts of 5v± (dashed lines) and 5B±/d4:iip (solid lines), at bottom left the upgoing and 
downgoing wave energy fluxes, at top right the profile of the Alfven speed, and at bottom 
right the profiles of OVa/Oz (dashed lines) and a/|Va<9 2 Va/<9z 2 | (dash-dot line), again all for 
a wave of angular frequency 0.07 rad s _1 (horizontal dotted line in the bottom right panel) 
on a coronal loop with length 75,000 km and magnetic field 15 G. Expressions for cut-off 
frequencies, at which the Alfven wave becomes evanescent, differ among various authors, 
but are of similar magnitude to quantities like OVa/Oz and a/| V^<9 2 Va/<9-2 2 |. Also plotted 

as a solid line is d (JWd dz) 2 /4 + \VAd 2 VA / dz 2 \ / 2, which is given as the cutoff frequency 



by iMusielak et al.l (119921 ) and iMoore et al.l (119911 ). As can be seen, most wave reflection 



occuring at the steep chromospheric density gradient between 2100 and 2200 km altitude 
will not be greatly affected by changing of the Alfven wave frequency in the range 0.01 - 0.1 
rad s _1 . Also, an increased wave frequency due to increased coronal magnetic field will not 
change its reflection properties, because the critical frequency also increases with Va- 

So while the energy transmission depends primarily on the loop length, and hence 
wavelength of the Alfven wave, from Figure 4, it appears that the magnetic field is more 
important in controlling the degree of fractionation. Weaker magnetic field means that the 
altitude where the plasma /3 ~ 1.2 is higher, and the depth of chromosphere beneath a loop 
footpoint where fractionation may take place is reduced. The slow mode wave generation 
is also stronger in weaker magnetic fields, and so as the magnetic field becomes weaker, the 
more ponderomotive acceleration can be saturated by the associated longitudinal pressure. 
Additionally, compared to waves precisely on resonance, waves slightly off-resonance give 
a very similar ponderomotive acceleration profile, but generate more slow mode waves in 
lower regions of the chromosphere, again biasing fractionation further towards the top of 
the chromosphere. Thus we may expect the minimum He abundance to be just to higher 
frequency than the actual loop resonant frequency, as observed in our calculations. Further 
from resonance the spike in the ponderomotive force generates both positive and negative 
accelerations counteracting one another, such that essentially all fractionation occurs low 
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down in the chromosphere, and He is thus not strongly depleted. 

The importance in this example of the relative suppression of fractionation low in the 
chromospher e sugges t s that elements such as C or S might "echo" this behavior of He. As 



explained in iLamingj ( 120121 ). these elements can be fractionated in the low chromosphere, 
assuming the ponderomotive acceleration is present here, because H is dominantly neutral. 
This means that elements fractionate roughly in proportion to their ion fraction, and both 
C and S are approximately 90% ionized. Once H is ionized, fractionation is suppressed 
unless the charge state fraction is much closer to 1; for example Fe and Mg have ionization 
fractions of 0.9995 and 0.998 respectively. Thus when the ponderomotive acceleration is 
concentrated at the top of the chromosphere, and when He is strongly depleted, S and 
C should have relatively little fractionation, and correspondingly when the acceleration is 
present in the low chromosphere and He/O is relatively unaffected, S and C should display 
some fractionation. 

Figures 9 and 10 show C/O and S/O as a function of frequency. Both show the expected 
qualitative behavior with frequency and magnetic field, although the minimum in these ratios 
appears always to be at the loop resonance, as a steeply peaked dip. A small distance above 
resonance the fractionation is a function of magnetic field and loop length with the highest 
fractionations for the highest magnetic field and the shorter loop lengths. The fractionation 
peaks at some value and then declines far away from resonance. 



4. Variation of the FIP Effect with Solar Activity? 

4.1. The Helium Abundances 

Aellig et al. (2001) and Kasper et al. (2007) both found that the solar wind Helium 
abundance varied with wind speed and the solar cycle. Although there is a large scatter in 
the individual data points (Aellig et al. 2001) the He/H ratio averaged over a velocity bin 
increases almost monotonically with wind speed. The variation with wind speed is wider 
during solar minimum than during solar maximum in the sense that the strongest average 
He depletion is seen at solar minimum. 

We calculate fractionations relative to O, on the assumption that all elements can be 
treated as minor ions (this is not true for H). To be certain that this did not compromise our 
results we looked at the publicly available level 2 ACE SWICS He/O and He velocity data 
over the life of the mission so far. Figure 11 reproduces the findings of Aellig et al. (2001) 
and Kasper et al. (2007) with respect to wind speed, in that lower He velocities correspond 
to lower average He/O ratios in a mostly monotonic fashion. We do not find the strong 
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dependence of the helium abundance on phase of the sol ar cy cle. It must be rem embered 
that we are considering He/O, whereas lAellig et al.l ( 1200 ll ) and Kasper et al.l ( 120071 ) worked 
with He/H. 

These authors followed several antecedents and discussed their results in terms of inef- 
ficient Coulomb drag or gravitational settling. We are going to argue that this abundance 
anomaly is part of the same process that gives rise to the FIP effect, but before doing so 
we briefly review thi s alternative. Gravit ational settling is known to occur on the Sun in 
quiescent strea mers (IFeldman et al.l 119981 ) . where Fe is observ ed to be depleted relative to 



other elements. iNoci et al.l (ll997[ ). following iGeiss et all ( Il970l ). make the point that if such 
processes were the origin of the solar wind He abundance depl etion, oth e r hea vy ions should 
also be depleted, which is now known not to be the case. iBochslerl (120071 ) revisits this, 



comparing He, O, and Ne for which He appears most affected by gravitational settling, or 
inefficient Coulomb drag, in the coron a. In this cas e the He abundance should correlate with 
the proton flux, which is not observed (jWangil2008l ). He and other m inor ions are also known 



to flow faster than H in the solar wind obse rved in situ at 1 AU (e.g. iNeugebauer et al.lll996 



Kasper et al.ll2008t iBourouaine et al.ll201ll ) which conflicts with the notion of a proton driven 
wind dragging the rest of the minor ions out with it, unles s a second stage of wind acce lera- 
tion is invoked further out in the heliosphere (jWangl 120081 ) . lLaming fe Feldmanl ( 120031 ) also 
measure a depleted He abundance in coronal hole and quiet solar region at altitudes below 
that where He is presumed to be depleted in IBochslerl (120071 ). 



We believe we can better understand the variation of the solar wind He abundance 
in the context of our model of the FIP effect. Figure 4 shows the variation of He/O as a 
function of the frequency relative to the resonant frequency. The strongest He depletions 
are only seen just above the resonant frequency. He/O approaches its photospheric value as 
one moves further away from resonance. Consider a loop being heated by a process that also 
generates Alfven waves. So long as the perturbation to the magnetic field is small compared 
to the loop initial magnetic field, only the resonant frequencies of the loop would be excited. 
However the heating would evaporate material from the footpoints of the loop up into the 
loop, changing the resonant frequency. Low levels of activity and heating would only change 
the resonance by a small amount. More active regions would have greater heating, greater 
evaporation and thus move the loop further away from the original resonance. The He 
depletion due to the originally excited waves then depends on how far from resonance the 
loop became. Thus according to our parameter search, higher slow wind speeds should be 
associated with plasma that was originally in (shorter) loops with higher magnetic field. 
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4.2. The Sulfur and Carbon Abundances 

In the above picture of loop excitation, we always expect the excited wave to be slightly 
off resonance due to chromospheric evaporation. Thus from figures 9 and 10 we expect the S 
and C fractionations to be a function primarily of the magnetic field in the loop. This means 
we predict that one should see the average C and S varying with solar cycle and wind speed, 
just as is seen for He. While S is not available in the level 2 ACE products we can test this 



hypo thesis for C. In figure 12, we plot the average C/O relative to solar (0.550 lAsplund et al. 



20091 ) for 9 velocity bins as a function of time, in the same way as figure 11 for He/O. For at 
least the years 2000 to 2004 and during 2005 a variation with wind speed is apparent, with 
the slowest wind speeds showing a lower C/O fractionation. The period around 2003 also 
corresponds to a period of strong He/O fractionation with wind speed. This fits well with 
the prediction of our model. Further confirmation could be found if S/ O could be analy zed 



sec 



as well. A similar correlation between He/O and C/O has been found (IKo et al.l 12006 
their Figure 11), who analyzed ACE Level 2 data collected over a week timescale as fast 
wind from an equatorial coronal hole transitions to slow wind from the surrounding corona. 



4.3. The Coronal Origin of Alfven Waves 



Our finding that Alfven waves on resonance with a coronal loop give FIP fractionations 
in better agreement with observations than other frequencies reinforces our suggestion of a 
coronal origin for the waves. No chromospheric origin would naturally achieve such selec- 
tivity, unlike the expectation from a coronal origin. For example, consider a flux tube with 
cross-B density gradient scale length /, undergoing a kink mode oscillation which become 
Alfvenic at a resonant surface where the w ave frequency oj = 2L/Va, where L is the loop 
length and Va is the Alfven speed (see e.g. iRuderman fc Roberts! |2002| ). The width of this 
resonant layer is 5 ~ (lu/tu) 1 ^ 3 , where v = clVa/R is the kinematic viscosity, a is the loop 
radius and R is the Reynolds number. The Alfvenic velocity fluctuations are larger than 
those of the kink mode by a factor 1/5 = (PR/aL) 1 ~ 0.1i? 1//3 under typical conditions, 
and are constrained to be resonant with the flux tube and their location of excitation. 

The nanoflare paradigm suggests that the dominant loop footpoint motions are of much 
lower frequency, and do not excite oscillations in the loop but act so as to build up magnetic 
stresses in the corona. These stresses periodically release themselves, in what has become 
known as a "nanoflare", as a current sheet develops. Rappazzo et al. (2007, 2008) discuss 
the buildup of magnetic stresses within the framework of turbulence phenomenology, where 
it appe ars that velocity pertu rbations similar to the ~ 30 km s _1 observed should be ex- 
pected. lLongcope et al.l ( 120091 ) conjecture that in impulsive reconnection in post-flare loops, 
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only about 10% of the liberated magnetic energy is converted directly into heat, the rest 
reappearing as kinetic energy that ultimately drives turbulence. So long as the perturbation 
to the magnetic field is sufficiently smaller than the ambient magnetic field, the boundary 
conditions imposed by the loop footpoints will restrict the modes excited by the reconnection 
event to those that are resonant with the loop. 



Sturrockl (119991 ) gives a pedagogic review of the mechanisms by which various wave 



modes may be excited by reconnection. The reconnected field line is generally distorted, and 
this can either propagate away from the reconnection site as an Alfven wave, or emit mag 



netoacoustic waves traveling perpendicularly to the magnetic field direction. Kigure et al. 



( 120101 ) explicitly consider the generation of Alfven waves by magnetic reconnection, and find 
that a significant fraction of the magnetic energy released (several tens of %, depending 
on geometry and plasma /3) can be carried off by Alfv en or mag n etoac oustic (fast or slow 
mode) waves, with Alfven waves dominating for < 1. iLiu et al.l (120111 ) discuss the role of 
temperature anisotropics and wave generation by the firehose instability in the outflow. 



5. Origin of Slow Speed Solar Wind 



Based on the variation on the helium abundance in the slow speed wind wind, we argue 
that the faster slow speed wind should have originated in magnetic flux tubes with higher 
magnetic field, with the slower components of the slow wind arising from flux tubes with 
weaker ma gnetic field. Slow wind s peed is known to correl ate with the expansion rate of the 
open field (jWang fc Sheeleyll 19901 ) . IWang fc Sheeleyl ( 120031 ) argue that at sunspot maximum, 
more slower speed slow wind is produced, and that this comes mainly from small sheared 
open field regions located near active regions. These have relatively high footpoint field 
strengths, and correspondingly high expansion factors. The high coronal heating rate near 
the footpoint (heating is presumed to correlate with magnetic field strength) increases the 
solar wind density, thereby reducing its speed f or the sa me energy input . Our inference 
concerns not this open field considered by IWang (120081 ) and iBochslerl (120071 ) , but the closed 
field with which it presumably reconnects releasing plasma into the solar wind. 

Attempts have been made to map in situ slow solar wind observations by Ulysses and 
ACE spacecraft back to a magnetic field sou rce surface at a heliocentric distance of arou nd 2.5 
R Q , and thence back to the solar disk itself ( jNeugebauer et al.ll2002t iLiewer et al.ll2004j ). with 
the result that at sunspot maxim um (bu t not a t activity minimum) , open field lines in active 
regions are the most likely source. IWangi (120081 ) studied the solar wind helium abundance and 
found that it correlated best with the source-surface magnetic field strength, which was found 
to be strong in fast wind and slow wind at sunspot maximum, but weak in the low-speed 
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solar wind emanating from coronal hole boundaries at solar minimum. Increased proton flux 
in the low-speed solar wind should i ncrease the H e abundance if inefficient Coulomb drag is 
important, but as mentioned above, IWang J2008h found essentially no correlation between 
these variables. 

If the magnetic field in the originally closed loop correlates with the strength of neigh- 
boring open field with which it must reconnect, our conclusions match these results. From 
Figure 3, the dependence of He/O on magnetic field for a given loop length appears to be 
much stronger than the dependence on loop length for a given magnetic field, although there 
is a tendency for longer loops to have stronger helium depletions. 

Our results strongly suggest that the solar wind helium abundance is set in the chromo- 
sphere , along with the other abundance modifications that result in the FIP effect. iByhring 
(120111 ) reaches a similar conclusion, but argues that it must arise from gravitational settling. 
However such processes; diffusion, gravitational settling and inefficient Coulomb drag are 
too slow to have an appreciable effect in a solar atmosphere that is increasingly being ap- 
preciated as a dynamic environment . For example time scales in diffusion models ra nge from 
tens of hours (IHansteen et al.l 1 1994 ) to days or weeks (IKillie fc Lie-Svendsenl 120071 ) , and so 
such possibilities must be considered highly implausible. 

Our calc ulations necessarily have approximations in the treatment of wave physics (see 
Lamingll2012l . for fuller discussion). However all of the phenomena we discuss have qualitative 
explanations in terms of the structure of the chromosphere and its ionization balance. While 
absolute numbers for fractionations may vary slightly in future works as chromospheric 
models and atomic physics improve, we expect that the various correlations between element 
abundances and the relation of the FIP fractionation to the wave frequency in terms of the 
loop resonant frequency to be robust. Most of the variations we have discussed center 
around the "spike" in the ponderomotive acceleration at the top of the chromosphere where 
the density gradient is steep. At first sight it might seem a coincidence that this "spike" 
occurs where the hydrogen in the chromosphere is becoming ionized, and thus allows for 
FIP fractionation. This is not the case. The steep density gradient arises precisely because 
hydrogen is becoming ionized, and so cooling by radiation in Lyman a is becoming quenched. 
Consequently heat input into the plasma cannot be radiated away, and the temperature must 
rise. The density decreases to maintain constant pressure, and waves in this region give rise 
to the strong ponderomotive acceleration. 
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corona 




chromosphere 



Fig. 1. — Cartoon illustrating the model. The integration of the Alfven wave transport 
equations is begun in chromosphere A with one downgoing wave. Integrating back through 
the corona to chromosphere B, forward and backward propagating waves develop. At chro- 
mosphere B, the wave is either reflected back into the corona, or transmitted into the chro- 
mosphere, depending on the wave frequency with respect to the resonant frequency of the 
loop. Resonant waves reflect, nonresonant waves precipitate down. 
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Fig. 2. — Depiction of force-free chromospheric magnetic field for coronal loop model with 
B = 15G, calculated from Athay (1981). Solid lines indicate magnetic field lines. Dashed line 
indicate contours of the Alfven speed. The plasma f3 — 1.2 layer is at 700 km altitude above 
the photosphere, and moves up and down for weaker or stronger coronal fields respectively. 
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Fig. 3. — Fractionations relative to oxygen relative to their photospheric values as a function 
of FIP for the case of 15 G and 75,000 km. Overlaid are the measured fractionations of the 
elements included in von Steiger et al. 2000. The strongest He depletion corresponds to the 
purple line, illustrating the resonant case. As one moves off resonance, the He abundance 
increases. Other elements are predicted to vary either in correlation (e.g. Al, Mg, S, C), or 
anti-correlation (e.g. Rb, W) with He. 
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Fig. 4. — He/0 relative to their solar ratio as a function of the frequency divided by the 
resonance frequency. For 5G and 100,000 km the points where it was not possible to obtain 
Fe/O of 4 for any amplitude, and these point are omitted. For 10 G and 75,000 or 100,000 
km, it was also not possible to achieve Fe/O ~ 4 for any amplitude, and these points are 
also omitted. A trend of increasing He depletion with decreasing magnetic field and possibly 
also increasing loop length can be seen. 
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Fig. 5. — Amplitude in km s _1 at the peak of the loop as a function of frequency relative 
to the resonance frequency. The absent points have the same meaning as in Figure 3. All 
curves are selected to give Fe/O close to 4, but as can be seen, close to resonance, this results 
in very similar wave amplitudes, even though loop lengths and background magnetic fields 
are varying. 
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Fig. 6. — The energy transmission coefficient as a function of wave frequency (expressed in 
terms of the loop resonant frequency) for the loop lengths and magnetic fields considered in 
this paper. A transmission coefficient of unity is achieved on resonance for all loops. Away 
from resonance, shorter loops (with higher frequency waves) transmit more energy between 
the chromosphere and the corona. The magnetic field variation is less significant. 
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Fig. 7. — Chromospheric section of the model loop, length 75,000 km, magnetic field 15 
G, angular frequency 0.07 rad s _1 . The upper left panel gives the density (solid line) and 
temperature (dashed line). Bottom left gives the ponderomotive force and slow mode wave 
(i.e. longitudinal wave excited by the Alfven waves) amplitude. Upper right; ionization 
fractions of Fe, Mg, S and C. Bottom right; fractionations of Fe, Mg, S, and He, as well as 
the ionization fractions in grey of H (dash), C (dash triple dot), and He (dash dot). 
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Fig. 8. — Chromospheric section of the model loop, length 75,000 km, magnetic field 15 
G, angular frequency 0.07 rad s -1 . The upper left panel gives the amplitudes of Elsasser 
variables in km s _1 (5B /y/Anp solid lines, 8v dashed lines, real parts in black and imaginary 
parts in gray.) At lower left are shown the wave energy fluxes in erg cm _2 s _1 for left and 
right going waves. The thin dotted line shows the difference in energy fluxes divided by the 
magnetic field strength and should be a horizontal line if energy is conserved. Upper right 
gives the Alfven speed Va profile through the chromosphere. The bottom right right gives 
the variation of OVa/Oz (dashed line) and ^VaO^a/Oz 2 (dot-dashed line). The solid line 

gives J(dV A /dzy /4 + \V A d 2 V A /dz 2 \ /2, which is the Alfven wave cutoff frequency as given 
by Moore et al. (1991) and Musielak et al. (1992). 
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Fig. 9. — S/0 relative to their solar ratio as a function of the frequency divided by the 
resonance frequency, similar to Figure 4. 
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Fig. 10. — C/O relative to their solar ratio as a function of the frequency divided by the 
resonance frequency, similar to Figure 4. 
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Fig. 11. — ACE SWICS yearly average of the measured 2 hr averaged He/0 ratio relative to 
photospheric (126) from 1998 to 2010 as grouped by He velocity in 50 km s _1 bins. Errors 
are the standard error on the mean. A trend of higher He/O abundance ratio with higher 
solar wind speed is clearly visible, though strong variation with phase of the sunspot cycle 
is not seen. 
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Fig. 12. — Same as Figure 11, but for C/0 relative to photospheric (0.550). A similar 
correlation of C/0 with wind speed is seen, especially around solar minimum years 2002 - 
2004, and 2006. 



